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We review recent analytic and numerical results concerning the confinement scenario in Coulomb 
gauge. We then consider a local, renormalizable, BRST-invariant action for QCD in Coulomb 
gauge that contains auxiliary bose and fermi ghost fields and sources. When the auxiliary fields are 
integrated out, one obtains the standard Coulomb gauge action with a cut-off at the Gribov horizon. 
We use the local formulation to calculate the leading correction to the Stefan-Boltzmann equation 
of state at high temperature due to the cut-off at the Gribov horizon. It is of order g 6 , which is 
precisely the order at which the infrared divergence found by Linde divergence first occurs. No such 
divergence arises in the present calculation because the propagator of would-be physical gluons is 
suppressed in the infrared due to the proximity of the Gribov horizon in infrared directions. 

PACS numbers: 12.38.Aw, 12.38.Mh, 12.38.-t, ll.10.Wx, 11.15.-q, ll.10.Gh 



I. INTRODUCTION 



It has poles at 



A. Confinement scenario in Coulomb gauge 

One would like to explain the presence of a long-range 
force that confines colored objects. At the same time 
one is faced by the apparently contradictory requirement 
that the massless gluons that are supposed to transmit 
this force are absent from the physical spectrum. 

This paradox is addressed in a scenario in Coulomb 
gauge that was originally developed by Gribov [l|. In 
this scenario, the 3-dimensionally transverse, would-be 
physical gluon propagator, DAiAj (k, k^), is suppressed at 
small k because of the proximity of the Gribov horizon 
in infrared directions, so physical gluons are absent from 
the physical spectrum. On the other hand confinement 
of quarks, or any colored object, is explained by the long 
range of the instantaneous part, V cou \{R), of the time- 
time component of the gluon propagator, 



g 2 D AoAo (x,t) = Koul(|x|)5(f) + P(x,t), 



(1) 



fc 2 = £ 2 ( k)^k 2 + ^ 



(3) 



in agreement with the Coulomb-gauge energy obtained 
by Gribov. 

All states are physical in Coulomb gauge, and the 
instantaneous part, V cou \(R), of the gluon propagator 
is directly related to the energy of a quark pair state, 
\^qq) — g(0)g(J?) I^o)- This state in Coulomb gauge con- 
tains massive quark-anti-quark sources at separation i?, 
and is the vacuum state of pure glue. tyq q is invariant 
under the remnant gauge symmetry of time-dependent, 
space-independent gauge transformations g(t) that is left 
unfixed by the Coulomb gauge condition, 



^cU,(x,i) = 0. 



(4) 



Let L(x,T) be the time-like Wilson (not Polyakov) line 1 



that couples universally to color charge. Here P(x, t) 
is a non-instantaneous vacuum polarization term that 
is screening, whereas the color-Coulomb potential, 
V cou i(R), is anti-screening 0]. The concentration of prob- 
ability at the Gribov horizon causes V^ ou i(i?) to have the 
long-range, confining property lim^^oo V^ou^i?) = oo. 
The different behavior of space and time components 
of the gluon propagator is possible because manifest 
Lorentz invariance does not hold in Coulomb gauge. 

The local, renormalizable theory which will be pre- 
sented here yields results in agreement with this scenario. 
For example, at tree level, the 3-dimensionally transverse 
gluon propagator, 



D(k, ho) 



(fc 2 + k 2 )k 2 



(2) 



L(x, T) = P exp / dx Q gt b A b (x, x ) 



(5) 



for quarks in an irreducible representation r of SU(N), 
of dimension d r . Here t a is a basis of the Lie algebra 
of SU(N), [t a ,t b ] = f abc t c in the representation r. The 
correlator of two Wilson lines, 



G(R,T) = d^(Trp(x,T)L(y,T)}, 



(6) 



may be expressed in terms of the Hamiltonian H and the 
quark-pair state, 



G(R,T) = (* qq \ e -("- £ o)T 



(7) 



is strongly suppressed in the infrared, vanishing like k 2 . 



1 Here we follow the discussion of 
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where Eq is the vacuum energy, and R = |x — y|. We 
have 



G(R,T)=Y. C " 



(8) 



sees from © that, for large T, e(R,T) converges to the 
minimum energy of a quark pair in the representation r 



lim e(R,T) = V r (R) + E s 



(15) 



where e n — E n — Eq is an excitation energy and c n 



K*»l*5,>f 

derivative, 



is a positive coefficient. The logarithmic 



e(R,T) 



d\n[G(R,T)] 
~~ df 



(9) 



yields the mean excitation energy above the vacuum en- 
ergy of the quark-pair state e~( H ~ E ° %>T t 2 '$> q q . 
For small T, the Wilson line is given by 



L(x,T) = l+.g J dx t b A b {x,x Q ) 



(10) 



and 



G(R,T) = l-TE se + (g 2 /d r )Tr(t b b c ) (11) 

T T 

x J J dx dyo{A b Q (x,x )A^(y,yo)) + ... , 



where E se is the regularized self-energy of the two quarks. 
It might be thought that the integral is of order T 2 . How- 
ever, with 

(A b (x,x )AZ(y,y )) = 6 bc D AoAo (x-y,x -y ), (12) 

the instantaneous term in fT]). proportional to 8(xq — yo), 
gives a contribution of order T, and we obtain, for small T 



G(R, T) = 1 - T[E se - C r V coul (R)} 



(13) 



where G r = —t b t b > is the Casimir invariant in the 
representation r of SU(N). This gives, 



lim e(R,T) = E se - C r V coul (R). 



(14) 



[The minus sign in front of V cou \(R) occurs because, 
by (JTJ) and ([12]), V cou \(R) is the color-Coulomb poten- 
tial between like charges.] Thus the physical energy of 
the quark-pair state ^ qq is given by the instantaneous 
part of the gluon propagator in Coulomb gauge. Con- 
sistent with this, the color-Coulomb potential V cou \(R) 
and, more generally, gAg are renormalization-group in- 
variants 

It may be shown that the remnant gauge invariance 
of time-dependent gauge transformations g(t) is sponta- 
neously broken, or not, according as V cou \(R) is not, or 
is, confining 

With these results, it is simple to establish that there is 
"no confinement without Coulomb confinement" p| ■ One 



which, for a confining representation, is the energy of the 
flux-tube ground state plus E se . One also sees from © 
that the time derivative of e(R, T) gives the negative of 
the variance of the excitation energy 



8e(R,T) _ 



dT 



(c 2 



2 ) <o, 



(16) 



which is negative. Thus e(R, T) is a monotonically 
decreasing function of T, and so e(R, 0) > e(R,oo). 
We conclude from (fT4)) and (p~5|) that if the quark- 
pair energy V r {R) in some representation r is confining, 
lmifj^oo V r (R) — * oo, then 



V r (R) < -C r V coul (R) 



R 



(17) 



Thus if the quark-pair energy is confining for some rep- 
resentation r, then also the color-Coulomb potential is 
confining —V cou \(R) — > oo. This analytic result strongly 
supports the confinement scenario in Coulomb gauge. 



B. Recent numerical results 

Numerical studies provide a valuable laboratory for 
testing various confinement scenarios. In accord with 
the above scenario, it was found that the equal-time 
would-be physical gluon propagator Da^j (k)|t = o is in- 
deed suppressed at small k, while the fourier transform 
K ou i(k) is enhanced at small k, which corresponds to 
a long-range, color-Coulomb potential V CO u\(R). It was 
subsequently found numerically 0] that the long-distance 
behavior of the color-Coulomb potential is consistent 
with a linear increase at large R, 



rs " 1 ^coul 

R. 



(18) 



with the Coulomb string tension given approximately by 
<y coui ~ 3er. Here a is the physical string tension between 
a pair of external quarks as determined from a large Wil- 
son loop in the fundamental representation. Other stud- 
ies provided additional support for the confinement sce- 
nario in Coulomb gauge, including in particular the con- 
sistency of Coulomb-gauge and center-vortex scenarios 
[H, Q . These numerical studies were certainly encourag- 
ing for the confinement scenario in Coulomb gauge. 

However, as in the old tale of the sorcerer's appren- 
tice, it was unexpectedly discovered that this "confine- 
ment scenario" works also in the deconftned phase. In- 
deed it was found numerically that the long-distance be- 
havior of V^ou^i?) is consistent with a linear increase, 
Ccoui > 0> at temperatures T above the phase transi- 
tion, T > T Cl where the physical string tension vanishes, 
(7 = [3j. Investigation of the temperature dependence 
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of er cou i revealed that at high T in the deconfined phase, 
the Coulomb string tension increases with T, consistent 
with a magnetic mass [9|, 

<^uiCn -cg 2 (T)T. (19) 

Thus, from the numerical evidence, what has been called 
the Gribov confinement scenario works as well in the de- 
confined phase of QCD as in the confined phase. 

Although this came as a surprise, it was a surprise 
that could have been predicted. Recall that the tem- 
perature T determines the extent, f3 = T -1 , of Eu- 
clidean space-time in the time direction. On the other 
hand gauge fixing to the minimal Coulomb gauge, de- 
fined in (PiTj) below, is done at each t, so the Coulomb 
gauge condition (j4|) holds at each Euclidean time t and 
the cut-off at the Gribov horizon applies to 3-dimensional 
configurations A(x, t) at each time, t, regardless of the 
extent (3 in the time direction. Thus the arguments which 
lead to suppression at small k of the space components 
DAfAjO^, ko) of the gluon propagator and enhancement 
at small k of the time component Da ,a {^ fco) hold at 
all temperatures T, including in the deconfined phase. 

It thus appears, from both numerical and analytic evi- 
dence, that what originated as a scenario for confinement 
provides a more general framework for QCD that holds 
both in the confined and deconfined phases. Although 
this is surprising at first, there is no paradox because the 
bound (|17[) provides a necessary condition for confine- 
ment, but not a sufficient condition. 

(To see that the converse may be false — namely that 
the color-Coulomb energy may diverge even though the 
quark-pair energy is finite — suppose that the ground- 
state quark-pair energy is confining when the quarks are 
in the fundamental representation, V{ nn d(R) ~ R. In this 
case according to the bound ([T7j). —V CO ui(R) oo. But 
then the color-Coulomb quark-pair energy — C r V cou i(R) 
is confining in every representation r with C r ^ 0. This 
includes the adjoint representation which however is not 
confining.) 

C. Present results 

In a first attempt to account for the influence of the 
cut-off at the Gribov horizon on the equation of state 
of the gluon plasma, a simple Planck formula was taken 
as Ansatz ([Hi), with a dispersion relation given by the 
Gribov energy E = (k 2 + ^p-) 1 / 2 , where the Gribov 
mass m was assumed to be temperature- independent. 
In the present work we develop a systematic calcula- 
tional scheme in which the cut-off at the Gribov hori- 
zon is effected within the framework of a local, BRST- 
invariant, renormalizable quantum field theory, and a 
renormalizable gap equation. The equation of state is no 
longer given by a simple Planck formula, and the Gribov 
mass becomes a calculable quantity that is temperature- 
dependent, m = m(T). We calculate the leading correc- 



tion to the Stefan-Boltzmann equation of state at high 
temperature due to the cut-off at the Gribov horizon. 

In 1980 Linde showed that standard and resummed 
finite-temperature perturbation theory suffer from in- 
frared divergences [11(. Since then no solution has been 
found, although the infrared divergences may be avoided 
by introducing a magnetic mass for the gluon in an 
ad hoc manner. Some time ago it was suggested that 
these divergences arise because standard perturbation 
theory neglects the suppression of infrared gluons re- 
quired by t he p roximity of the Gribov horizon in infrared 
directions [12j |. No infrared divergence appears in the 
present calculation because the tree-level propagator @ 
of would-be physical gluons is strongly suppressed in the 
infrared. The correction obtained here is of order g 6 , 
which is precisely the order at which the Linde diver- 
gence would otherwise occur. 

We show by integrating out auxiliary fields that the lo- 
cal, renormalizable formulation presented here is equiv- 
alent to the non-local action and gap equation that had 
been obtained previously from the standard Coulomb- 
gauge action by imposing a cut-off at the Gribov hori- 
zon (l3| . Determination of the number of independent 
renormalization constants (there are two) will be pre- 
sented elsewhere [i~4| . 



D. Other approaches 

Equal-time correlators in Coulomb gauge may be cal- 
culated using the Schwinger-Dyson e quat ions of the 
Hamiltonian formulation [19, 

mm mill- 

The gluon 

energy obtained by variational calculation accords at 
both high and low momentum with the Gribov energy 
©, provided that the Faddeev-Popov determinant is 
properly accounted for [T^, Oil]. Indeed the Faddev- 
Popov determinant dominates the low-energy dynamics. 
A systematic study of the infrared limits of propagators 
and vertices has been reported recently in [19j , which 
incorporates earlier results [2(| [Hj]. The Schwinger- 
Dyson equations in Coulomb and Landau gauge should 
be quite reliable in the infrared limit because the weight 
in the infrared limit is precisely given by the Faddeev- 
Popov determinant, det(— Didi) with negligible contri- 
bution from the gluon wave- functional or action (22j, a 
phenomenon known as ghost dominance [H, (24[. This 
is used to advantage with the horizon condition and 
non-renormalization of ghost-gluon vertices and, as a re- 
sult, the infrared limit of the Schwinger-Dyson equations 
in Coulomb and Landau gauges decouples from finite- 
momentum correlators. Moreover errors from trunca- 
tion of 3- vertices have been controlled [HI, [2f| By 
contrast the semi-perturbative (or semi-nonperturbative) 
method described here may not give accurate results in 
one-loop approximation for infrared quantities. Never- 
theless both methods have their advantages. In partic- 
ular the present approach also gives the correlators in 
Coulomb gauge at unequal times and finite temperature, 
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and does not rely on a variational Ansatz for the wave- 
functional that would be needed to go beyond the in- 
frared limit. 

The instantaneous color- Coulomb potential has re- 
cently been used in Dyson-Schwinger and Bethe-Salpeter 
equations to find solutions for pseudoscalar and vector 
mesons [13] ■ Cancellation of the energy divergences of 
Coulomb gauge has been demonstrated at the two-loop 
level [28|, and renormalization and cancellation of en- 
ergy divergences to all orders in perturbation theory in 
Coulomb gauge has been elucidated recently, [29l. [30j. 

Equivalence of a cut-off at the Gribov horizon to a 
modified local action with auxiliary fields was first es- 
tablished in Landau gauge some time ago [3l[. It was 
then shown that the horizon condition renormalizes con- 
sistently, and that the ghost propagator in Landau gauge 
has a ™ or dipole singularity 32j. The symmetries 
of the local theory in Landau gauge were exhibited and 
algebraic renormalizability established [HI, H3 ■ More re- 
cently, using this local action in Landau gauge, the gap 
equation was determined to two-loop order, and it was 
verified at the two-loop level that the gho st p ropaga- 
tor has a dipole singularity in the infrared [35|. It has 
also been found to one-loop order that the gluon prop- 
agator in Landau gauge vanishes like fc 2 , and that the 
renormalization-group invariant coupling a s (k ) ap propri- 
ate to the Landau gauge is finite at fc = [36| . Thus 
the elements of Gribov's scenario in Landau gauge [lj 
have been derived from a local, renormalizable, BRST- 
invariant action. The Landau gauge case has been re- 
viewed in [37| . 

The Coulomb gauge provides a more straightfor- 
ward confinement scenario than the Landau gauge (but 
see [11]), and for this reason Gribov and others turned 
to the Coulomb gauge. The Coulomb gauge is also 
well suited to finite-temperature calculations because the 
heat-bath provides a preferred Lorentz frame, so sym- 
metries of the Coulomb gauge, which breaks manifest 
Lorentz invariance, are the physical symmetries of the 
system. Unitarity is immediate in the non-perturbative 
Coulomb gauge because of the equivalence to a canon- 
ical system whereas it is problematical in the non- 
perturbative Euclidean Landau gauge. 



E. Organization of paper 

In sect. II we briefly review previous work on which 
this article relies, which is scattered in a number of dif- 
ferent places. A local, renormalizable, BRST-invariant 
action with auxiliary ghosts and sources is introduced 
in sect. III. In sect. IV the non-perturbative Coulomb 
gauge is defined by a gap equation that determines a non- 
zero "physical value" for sources. Unitarity and confining 
properties of the non-perturbative Coulomb gauge are es- 
tablished in sect. V. Free propagators are calculated in 
sect. VI. The gap equation is evaluated in one- loop ap- 
proximation in sect. VII and solved in sect. VIII. The free 



energy is calculated in sect. IX, and in sect. X the lead- 
ing high-temperature correction to the equation of state 
from the cut-off at the Gribov horizon is obtained. Our 
conclusions are presented in sect. XL The auxiliary fields 
are integrated out in Appendix A, and which results in 
the standard Coulomb gauge with a cut-off at the Gri- 
bov horizon. In Appendix B the partition function in 
Coulomb gauge with a cut-off at the Gribov horizon is 
expressed in canonical hamiltonian form. 



II. HOW TO CUT OFF A FUNCTIONAL 
INTEGRAL (OLD STUFF) 

In Euclidean quantum field theory the functional inte- 
gral 

Z = [ dA 5{d l A l ) det[A(A)d 4 ] exp{-S YM ), (20) 



must be cut off at the boundary dA of the "fundamental 
modular region" A, this being a region on the gauge- 
fixing surface that intersects each gauge orbit once and 
only once. The integrand is the familiar Faddeev-Popov 
weight Coulomb gauge. 

In Coulomb gauge we may take A to be the set of con- 
figurations Ai (x) each of which is the absolute minimum 
on its gauge orbit of the Hilbert norm, 



A = {4(x):||A||<||*A||}. 



(21) 



Here Ai(x) is a 3-dimensional configuration, 9 Ai — 
g~ 1 Aig + g~ 1 dtg is its gauge transform by a 3- 
dimensional local gauge transformation <?(x), and 
||A|| 2 = J d 3 x ^(x)! 2 is the Hilbert square norm. This 
gauge choice is done for each time t. 

At a local or global minimum, the functional FA,[g\ = 
|| 9 ^4|| 2 is stationary at g = 1, and the second variation of 
-Fa[<?] is a positive matrix (all its eigenvalues are positive). 
This translates into the Coulomb gauge condition diAi = 
0, and the positivity of the 3-dimcnsional Faddeev-Popov 
operator, 



Di(A)di > 0. 



(22) 



These two properties characterize a set known as the Gri- 
bov region f2, 



n = { A: diAi = 0, and - D l {A)d % > }, 



(23) 



which is the set of configurations each of which is a local 
minimum on its gauge orbit of the Hilbert norm. The set 
of local minima is larger than the set A of global minima, 

a c n. 

We do not possess an explicit description of the set A 
of global minima. However it has been argued from 
stochastic quantization [39j, a geometrically correct but 
unwieldy quantization of gauge theory, that integration 
over Q or over A gives the same expectation-value for 
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all ?7-point functions with finite n. This is consistent 
with abelian or center-vortex dominance scenarios, be- 
cause abelian and center-vortex configurations lie on the 
common boundary of the Gribov region and the funda- 
mental modular region Q. Either as exact truth, or to 
provide a model worthy of investigation, we replace the 
fundamental region A by the Gribov region f2 in the func- 
tional integral, 



Z 



dA S(diAi) det[-Di(A)di] exp(-S yM ). (24) 



The boundary <9f2 of the Gribov region, known as the 
"Gribov horizon" , is characterized in Euclidean Landau 
gauge or Coulomb gauge at fixed time (with D — > D — 1) 
by "horizon function" [311 ], 

G(A) ^fd^x [ gf abc A b (M- V ) ad gf dec A\ 

-(N 2 - 1)(D-1)}, (25) 

where M ab = -Df(A)d l is the D - 1 dimensional 
Faddeev-Popov operator, and M ah {M~ 1 ) hd f dec A\ = 
j aec A\. The horizon function vanishes on the boundary, 
G(A) = for A e dft, and is negative inside, G(A) < 
for A £ f2. Thus the partition function may be written 

dA 0[-G{A)] J(9 i A)det[-A(A)a i ]exp(-S'rM) I 

(26) 

where 6(—G) is the step function, and a product over 
all times t is understood. Configuration space with, for 
example, a lattice cut-off is a high dimensional space, and 
classical statistical mechanical arguments apply. Entropy 
favors population at the boundary <9f2 where G(A) = 0, 
and in the last integral we may make the replacement 
9[-G{A)] -► S[G(A)}. 2 The horizon function G(A) is 
a bulk quantity, like the hamiltonian H(A) in classical 
statistical physics. Just as the microcanonical ensemble 
is equivalent to the canonical ensemble, 5[H(A) — E] — » 
exp[— (3H(A)\, where [3 is determined by the condition 
(H) = E, likewise the microcanonical weight <5[G(A)] 
is equivalent to the Boltzmann weight exp[—jG(A)], so 
partition function is given by [l|, [l3|, [3l[ 



Z = j dA 5{d l A i ) det(-D z (A)d t ) exp(-S YM - jS h ), 

(27) 

where 7 is determined by the horizon condition 

(S h ) = 0. (28) 

Here, as a result of the product over all t, Sh is the inte- 
grated horizon function 



S h 



dt G[A(t)], 



(29) 



2 As an example of this, the radial probability distribution inside 
an ^-dimensional sphere of radius R is given by r N ~ 1 dr. This 



gets concentrated on the surface of the sphere, at 
N — > 00, thus f^dr r^- 1 -> J dr S(R - r). 



R, for 



where G(A) is given in (|2"5j) . Because (G[A(t)]} is inde- 
pendent of t, (|28|) is equivalent to (G[A(t)]) = 0. In terms 
of the free energy, W = In Z, the horizon condition may 
be written 



dW 

$7 



0. 



(30) 



The action Sh is non-local in space, but local in time. 
It may be expressed in terms of a local, renormalizable 
action by integration over auxiliary fields^as has been 
known for some time in the Landau gauge |31[ . Precisely 
the same argument holds for the Coulomb gauge. It is 
written out in Appendix A, where we start with the local, 
renormalizable action that is given below and integrate 
out the auxiliary fields to obtain l[27j) and (|30|1 . 



III. PERTURB ATIVE COULOMB GAUGE 

A. Faddeev-Popov action 

For simplicity we shall be interested in pure SU(N) 
gauge theory at temperature T, but there is no obstacle 
to extending our considerations to include quarks. It is 
described by a Euclidean action which, for pure SU(N) 
gauge theory, is of the form 



S = S Y m + sE 



(31) 



where Sym = / dP x Cym is the Yang-Mills action, with 



£>ym — — F, 



F MV = d^Ay - d v A ll + gAf, x A„ 



(32) 



(33) 



The s-exact term sS is defined below. Here g is the 
coupling constant, and we use the notation for the Lie 
bracket (A x B) a = f abc A b B c , where f abc are the fully 
anti-symmetric structure constants of the SU(N) group. 
The color indices a, b, c are taken in the adjoint repre- 
sentation, a = 1, A^ 2 — 1. We shall generally suppress 
the color index, and leave summation over it implicit. 
Configurations are periodic in Xq, 



A^Xi^o) = A^(xi,x + /?), 



(34) 



with period (3 = l/T, where T is the temperature. The 
integral over xq always extends over one cycle, J dxo = 

dxo. We are in D Euclidean dimensions. Lower case 
Latin indices take values, i = 1,2, ...,D — 1, while Greek 
indices take values fx = 0, 1, D — 1. 

In the BRST formulation there arc, in addition to A^, 
a pair of Faddeev-Popov ghost fields c and c and a La- 
grange multiplier field, b, on which the BRST operator 
acts according to 



sA 



sc. 



ib; 



sc = —(g/2)c x c 

sb = 0. (35) 
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It is nil-potent, s 2 = 0. Here is the gauge-covariant 
derivative in the adjoint representation, D^c = d^c + 
gAf, x c. 

The choice of 5 is the choice of gauge. Physics is in- 
dependent of 5, provided that it provides a well-defined 
calculational scheme. The standard Coulomb gauge is 
defined by the choice 5 CO ui = J d D x C cou \ 

■^coul — ^£coul — sdicA-i 

= idibAi — dicDiC. (36) 

The Lagrange-multiplier field b imposes the Coulomb 
gauge condition diA{ = 0. 



B. Auxiliary ghosts 

We first introduce the auxiliary ghosts and related 
sources in the perturbative Coulomb gauge where their 
role is trivial and superfluous. We shall then use them 
to write a local theory that is equivalent to the standard 
Coulomb gauge with a cut-off at the Gribov horizon. 

Recall that observables O are in the cohomology of s 
(namely s-invariant operators sO = 0, modulo s-exact 
operators, O ~ O + sX.) We may freely introduce addi- 
tional quartets of auxiliary ghost fields on which s acts 
trivially, because such fields cannot appear in the coho- 
mology of s, 



' B > 



B I 



SLU a B = 

sfe = 0. 



(37) 



Here a labels components in the adjoint representation 
of the global gauge group a = 1,...,N 2 — 1, and B is a 
mute index that is arbitrary for the moment. The fields 
6 B and (f>% are a pair of bose ghosts, while v B and u)% 
are fermi ghost and anti-ghost. The BRST method [4(| 
insures that physics is unchanged if we add to the ac- 
tion an s-exact term, S aux — J dPx C aux — J dPx s£ aux 
that depends on the auxiliary ghost fields. The auxiliary 
lagrangian, adapted to the Coulomb gauge, is taken to 
be 3 

Caux = S diU>%(Di<t>B) a 

= di4> a B {D i( j> B ) a (38) 
-d&%[ (Awb)° + (sAcx^) ]. 

The action and Lagrangian density are now 
S = J d D x C 

C = CyM + Aioul + Ajux) (39) 



where £ C oui and C aux are given in (|36|) and 

This action is renormalizable by power counting. Sym- 
metries and renormalizability of the analo gou s action in 
Landau gauge were established in [32| and [33[ , and sim- 
ilar considerations hold also in Coulomb gauge [lil ]. 



C. Auxiliary sources 



To derive the Slavnov- Taylor identities it is standard 
to introduce sources and L for the BRST-transforms 
that are non-linear (K Ul sA^) = (K^, D^c) and (L, sc) = 
(L,(—g/2)c x c). We also introduce auxiliary sources 
namely, pairs of bosonic sources, V„ a K an d V iB and pairs 
of fermionic sources N? B and N? B [32| . The extended 
action with all sources is defined by, 



d D x A 



A = C + K u sA u + Lsc- 



■DiwNi 



ViDi4> + sDiUbVi 
■NisDrf + ViVi-NiNi. (40) 

Here the color index a and the mute index B are summed 
over and suppressed, ViD^ = VisDi(j) B etc. We have 
also added the purely source term ViVi — NiNi, which 
will allow us to write the horizon condition in a multi- 
plicatively renormalizable form. 

D. quantum effective action 

The partition function is given in terms of the extended 
action by 

Z(J, K, L, V, V, N, N) — [ d$ exp[-E + (J,$)]. (41) 



where $ = (A,c,c,cp,ui,ui,ip) represents the set of all 
fields, and J the corresponding sources. The free energy 
is defined by W(J, K, L, V, V, N, N) = In Z, from which 
the quantum effective action T(<i>, K, L, V, V, N, N) is ob- 
tained by Legendre transformation. 

The quantum effective action satisfies the same 
Slavnov- Taylor identity in Coulomb gauge as in Landau 
gauge and the proof is the same as in Landau gauge , 

/ 6T ST 5T ST 
\~5KSA + SL 1c ~ 

56 SCo 



7 6T 

lb Tc 



N 



ST 



ST 



sv 



snJ 



(42) 



0. 



We will shortly assign non-zero "physical" values to the 
sources V and V . This gives an additional contribution 
—Vphj^i to the usual Slavnov- Taylor identity. 



More generally we may take the gauge-fixing term to 
be sf d D x d K c a K x an d the auxiliary action to be 

s J d x d K 

a n\ (D x rt> B ) a , where a KX is a positive symmetric 
matrix that may be diagonal 41] 



IV. NON-PERTURBATIVE COULOMB GAUGE 

In this section we describe a local theory which is 
shown in Appendix A to be equivalent to the standard 
Coulomb gauge with a cut-off at the Gribov horizon. 
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Modified local action 



B. Renormalizable horizon condition 



The mute index B on the auxiliary ghosts is now de- 
fined to be the pair of indices B = (c, ^t), where c labels 
components in the adjoint representation of the global 
gauge group, c = 1, N 2 — 1, and /i is a Lorentz index. 
With this assignment, the BRST operator acts according 
to 



SCO 

s 



ac - 0, 



(43) 



and the auxiliary lagrangian (|38p reads 

£ a «s = 0i<W„) a6 (44) 

-a^ fc [ (A^) Qb + ( 5 Acx^) ab ]. 

[The gauge-covariant derivative and the Lie commutator 
act on the first color index only, because the second color 
index is mute, thus (A0 M ) ac = d.^ + g(A, x M ) QC 
where {A l x c^) ab = f acd A\<f b \ We further stipulate 
that the "physical value" of the sources V and V is not 
zero but rather at 



f>a& Y7 a & „,l/2r rafa 



(45) 



where 7 is a constant with engineering dimension m that 
will be determined shortly. All other sources are set to 0, 
in particular N p h — N p \ l = 0. For purposes of a semi- 
perturbative expansion in g, described below, we write 



vV2 = 



(2iV)VV 



(46) 



and we take m to be of order g°. This yields the "physi- 
cal" local action 

S p h = ^\v=v ph ,v=v ph ,K=L=N=N=o = J d D x £ph (47) 



where 



YM 



^coul ~\~ £q 



(48) 



The specification of the non-perturbative Coulomb 
gauge is completed by requiring that m be a stationary 
point of the vacuum free energy 



ph 



or 



dm 



dm 



= 0. 



(51) 



This is a non-perturbative gap equation that deter- 
mines m as a function of Aqcd, and eventually also 
of the temperature T. As shown in Appendix A, this 
equation expresses the "horizon condition" obtained pre- 
viously [1 31 ] - It is compatible with multiplicative renor- 
malization because m renormalizes multiplicatively, as in 
Landau gauge [12]. 

From W p h = lnZ p h, and Z p h = J <i$ exp(— Sph), the 
horizon condition may be written 



dS, 



P h 



dm 



= 0, 



(52) 



where the expectation value is defined by the action S p h- 
By (gnj), this reads 



(D i (<l> i -^r-(gD i cXu )i r a ) = Z-^(N*-1)(D-1). 

(53) 

The second term vanishes, 



2 1 / 2 ™ 2 



(54) 



((gD iC x uj t ) aa ) = 0, 



because there is no cui term in the action (|39|) and (|48p . 
and translation invariance implies that the terms di<p and 
di4> do not contribute to (|53[) . The horizon condition 
reads 



2 1 /2 7 



(D-l)(iV 2 -l). (55) 



V. PROPERTIES OF THE 
NON-PERTURBATIVE COULOMB GAUGE 



A. Unitarity 



rn 



2Ng 



: (D-l)(N 2 -l) 



+ j2NW- 9 [D ^' S{D ^ )r - 



(49) 



and s(DiU)i) = + g(DiC x u)i). The vacuum free 

energy is the free energy W = InZ, with sources set to 
their physical value, 



Wp h (m) = -r ph (m) = -r| v=v 



,V=V„i 



(50) 



In Appendix B, the equivalence to a canonical system 
is established, and with it, unitarity. Unitarity is already 
manifest at tree level. Indeed the only poles in fco of 
the tree-level propagators (given below) occur at k 2 + 
k 2 + t^j- = 0, corresponding to the Gribov [l[ energy 

E = (k 2 + ^r) 1 / 2 . This real energy is consistent with a 
hermitian hamiltonian. By contrast, the tree-level gluon 
propagator of the non-perturbative Landau gauge [3l| 



has poles at k 2 



= 0, where k 2 = k 2 



This 



corresponds to complex energy E = (k 2 ± im 2 ) 1 / 2 . 
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B. A confining property 

Integration over Lagrange multiplier field b produces a 
functional J-function that imposes the gauge constraint 
diAi = 0. Suppose this integration is done, so Ai is 
identically transverse. The auxiliary field <p a ^ is also a 
Lagrange multiplier, and integration over it yields a func- 
tional (5-function that imposes the constraint 

- diDf<j>f - ^gf^A) = 0. (56) 

Quantization is done in a periodic box of finite spatial 
volume V . Integration of the last equation over V at fixed 
time t kills the first term, which is a spatial divergence, 
and yields a new constraint satisfied by Ai, 

J d 3 x A) = 0. (57) 
v 

This states that the zero-momentum component of Ai 
vanishes. This is an additional gauge condition on Ai 
that expresses a confining property: a gluon with 3- 
momentum k — is forbidden. In accordance with this 
condition, the propagator Da { a (k, fco), calculated be- 
low, vanishes with k. 



The leading term is of order g 2 , 

4 

S - 2 ^~iWgt (^-i)^- 1 )^ (6°) 

where the spatial quantization volume is V = L 3 , and 
the time extent f3 = T — 1 . Although this term is inde- 
pendent of the fields, it should not be ignored because, 
when the gap equation is solved for m = m(T), it gives a 
T-dependent contribution to the free energy. The terms 
in the action of order g , 

So = So,YAI + So,l + 5*0,2 + 5o, 3; (61) 

are quadratic in the fields and determine the "free" prop- 
agators, 

S ,ym = f d D x \{d„K - d v A%) 3 (62) 



5 ,i = / d D x (idib a A? - d i c a d l c a ) (63) 
So,2 = J d D x {d^fd^f - diQfd^f) (64) 



C. Automatic cut-off at Gribov horizon 

In Appendix A it is shown that after integrating out 
the auxiliary fields the action acquires a term propor- 
tional to the "horizon function" 

S h =Jd D x[ gf abc A b t (M~ 1 ) ad gf dec Al 

-(iV 2 -l)0D-l)]. (58) 

If one expands in eigenfunctions of the Faddeev-Popov 
operator M(A), the functional weight acquires the factor 
exp(— c 2 /A ), where Ao(A)is the lowest non-trivial eigen- 
value of M(A). By definition of the Gribov region fi, 
Ao (A) vanishes as A approaches the boundary <9f2 of the 
Gribov region from within. This constitutes an "auto- 
matic" , non-analytic cut-off of the functional integral at 
the Gribov horizon. 



VI. FREE PROPAGATORS 

We now develop a semi-pertubative calculational 
scheme. For this purpose we treat m as an independent 
parameter of order g°, and calculate all quantities per- 
turbatively, including the gap equation, to a given order 
in g. Then m = m(g, T) is determined by solving the the 
gap equation (|55[) non-perturbatively. 

We expand the action (jig]) in powers of g, 



So,3 = J d D x f ahc A\{^ ~ ^r. (65) 

The term 50,3 causes a mixing of the zero-order trans- 
verse gluon and bose-ghost propagators. 

To calculate the free propagators, we define the field 
that mixes with A\, 

^PJ^"^ (66) 

The orthogonal component ( 2 Af) 1/2 f abc { ( P^ a + Vf 1 ) an d 
other components of ip and Cp do not mix with Ai . More- 
over because of the Lagrange-multiplier term idjbAj, 
only the 3-dimensionally transverse part Aj of Aj con- 
tributes to propagators, and it mixes only with the trans- 
verse part ipj of tpj . Consequently, the free propagators 
of the fields Aj and tpj are determined by the mixed ac- 
tion 

S (A,ij) = Jd D x[ (l/2)[(ij) 2 + &AJ) 2 + [d^Jf] 
-im 2 Afipj ) , (67) 

which corresponds to the matrix in momentum space, 

/ k 2 + k 2 -im 2 \ 
y -im 2 k 2 J ' 

with determinant 



S — S-2 + S*o + 



(59) 



Ai = (k 2 + fc 2 )k 2 + m 4 . 



(68) 
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The free propagators are given by 



(2nY 

XPy(k)- 



(k 2 + fc 2 )k 2 + m 4 



(69) 



D Ai ^ (x-y) = 



(2tt) d 
xP„(k) 



T T^2exp[ik ■ (x - y)] 

u (70) 



(k 2 + fc 2 )k 2 



xP^-(k) 



D^(x-y)) = / T^TdTT T2^exp[ife- (a;-y)] 

k 2 + fcg 



(k 2 + fc 2 )k 2 



(71) 



Here fco = 2irn/f3 are the Matsubara frequencies, where 
n is any integer, T = 1/(3, and Pij(k) = Sij — kikj is 
the transverse projector. We have suppressed the trivial 
color factor 5 bc . The gap equation ([55)) reads 



i (A%0)^(0)) = — (D - l)(iV 2 - 1). (72) 



VII. GAP EQUATION IN ONE-LOOP 
APPROXIMATION 



When the left-hand side of the gap equation is eval- 
uated to zeroth order in g, using the mixed propaga- 
tor (1701). it reads 



d°- l k v-^ D-2 D-l 

r r E >u , ,2^.2 i „4 =H^r. (73) 



(2tt) z? - 1 ^ (k 2 + fc 2 )k 2 + m 4 Ng 2 



where we used Pa(k) = D — 2. 

To evaluate the sum over Matsubara frequencies, 



1 



feo 



(k 2 + fc 2 )k 2 +m 4 



J_ y L 

2 k 2 8 f-' 1 + (27m 



E 2 ^f3 n= ^ ±1 1 + {2^n/BEf 
where E = (k 2 + ^t) 1 ^ 2 , we use the identity 

sinh6» = 6»rj[l + (6i/n7r) 2 ] 



(74) 



(75) 



lnsinhfl = ln(9 + ^ ln[l + {8/mr) 2 }. (76) 



This gives upon differentiation 

cosh(9 _ 1 2 y, i 

sinh 9 ~ 9 + 9 1 + (mr/fl) 2 ' 

n— 1 v ' 7 

and we obtain for the sum over Matsubara frequencies 
1 cosh(8E/2) 



Q 



2k 2 E sinh((3E/2) 



1 



2k 2 E V exp(/3P) - 1 



(78) 



The first term is the value of Q at T = and the second 
is a Planck-type finite-temperature correction. The gap 
equation reads 



4 _ D f d D ~^k (D-2) 
* J (2TT) - 1 2k 2 E ' 



D - 1 



exp(/3P) - I J Ng 2 ' 

where we have made the substitution g 2 — > g 2 /i 4 
The first term on the left hand side, 



,2,,4-D 



4-D 



d^fc (D-2) 



(2ir) £) - 1 2k 2 P ' 



(80) 



converges for D < 4. The integral is readily evaluated by 
taking y — |k| 4 as integration variable, with the result 



I = 



m\ D -^D-2 1 r((D-2)/4)r((4-D)/4) 



2 (4n) D / 2 r((D - l)/2) 



In terms of e = (4 — D)/2, this gives 



(81) 



J = 



4tt 2 



i /m 2 \ , v i f ,, , r'(i/2)\ 
— in( — )+in(47r)+-(r , (i) + „,\ LM +i 



'(!)- 



r(i/2); 



(82) 

with neglect of terms that vanish with e. From the iden- 
tity HI, 



r(i/2) _ r(i) 
r(i/2) " T(iT 



2 In 2, 



(83) 



one obtains 
1 



/ = 



4tt 2 



1 / 777 2 \ 

- - In (— ) + ln(47r) - 7 - In 2 + lj , (84) 



where 7 = — T'(l) = 0.577215... . Upon making the 
standard MS subtraction, one obtains 



I 



4vr 2 ln ( 2 to 2 ) ' 



(85) 



After this substitution, the gap equation (I79|) reads 

1 3tt 2 



1 fefi 2 \ f dx 

4 in i2^j + y - 



exp(m/3u)-l Ng 2 (T) ' 



(86) 



71=1 



where u = (x 2 + -iy) 1 / 2 . 
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VIII. SOLUTION OF ONE-LOOP GAP 
EQUATION 

To solve the gap equation for to, we change unknown 
from to to to* = m/T so it reads 



IX. FREE ENERGY 



/(to*) = - m (—) + [— , 1 , 

2 \m* ) J u exp(TO*u) — 1 
o 



where 



y 



3?H 1, feu 2 (3 
In ' 



Ng 2 {T) 4 



y, (87) 



(88) 



To order g°, the free energy W = hiZ is given by 

expW = J d$exp(-S'_ 2 - So), (95) 

where S_ 2 , eq. (|60[) . is field-independent and of order 
g~ 2 , while So, eq. (j6"Tj) . is quadratic in the fields. We 
obtain 



W = W-2 + Wo 



(96) 



where, for D — 1 = 3 



The function /(to*) decreases monotonically, /'(to*) < 
0, with /(0) = oo and /(oo) = — oo, so this equation 
always has a unique solution 



W-2 



S-2 



3to 
2N? 



(N 2 - l)L 3 /3 (97) 



and 



m* = — = h(y). 



(89) 



expWo = J d$cxp(-S ). 



(98) 



Moreover by the inverse function theorem, h(y) is ana- 
lytic, so in this approximation there is no phase transi- 
tion. 

Asymptotically at high T, the running coupling g(T) 
in the MS-scheme is given by 



I 11 N . / 2ttT 



To evaluate Wo we observe first that the Faddeev- 
Popov ghost pair contributes a factor k 2 which is can- 
celled by the two factors of |k| _1 that come from the 
Ao and b integrations. Moreover all auxiliary bose and 
fermi ghost pairs with the same action give contribu- 
tions to W° that cancel. Each single 4-momentum mode 
of the action S°(A, ip), eq. (|6T|) . contributes A^ 1 ^ 2 to 



g 2 (T) 24 



— hi ( -■ ) , (90) expWo, where, by (J6SJ) , Ai = (k 2 + fc 2 )k 2 + to 4 . Cor- 

^Aw£j/ j. , . n _ _ n i . t_h j_i j_ _ • ^f) 



where A^gr is a physical QCD mass scale. We specialize 
to high temperature T where g{T) is small, so the one- 
loop gap equation should be a good approximation. We 
also suppose that at high temperature /i = /i(T) « 27rT, 



y(T) 



3tt 2 



Ng 2 (T) 



(91) 



According to (|87|) and (j88|) . large T or small g implies 
small to*, and for small to* we have 



dx 



1 



— / dx 



u exp(m*u) — 1 to* J x 



so the gap equation at high T simplifies to 



2 Km* 



2 3/2 



y(T). 



(92) 



(93) 



To leading order at high T or small to*, the first term 
on the left is negligible compared to the second, and we 
obtain, with to* = S and, with neglect of the second 
term in (pJTj) . 



responding to the bose ghost tpj that mixes with A b 
is an otherwise unpaired fermi ghost mode that con- 
tributes A^ 2 , where A 2 = k 2 . The net result is that 
for each 4-momentum mode of the A\ field we obtain 
a contribution (Ai/A 2 )~ 1/2 , where Ai/A 2 = A; 2 + E 2 , 

and E = (k 2 + ^r) 1 ^ 2 . There is an infinite product over 
all frequencies ko for each 3-momentum mode of the Ai 
field. For each 3-momentum k of the A\ field the result 
of this infinite product is the Planck partition function 
for a single mode 

oo 

exp[Wo(k)] = j^exp(-nE0) 

n=0 

= [l-cxp(-£/3)]-\ (99) 

Wo(k) = - ln[l - exp(-E0)]. (100) 

The sum over the A^ 2 — 1 color modes, the 2 degrees of 
transverse polarization, and over all 3-momentum modes 
k yields, with £ k -y j^, 

d 3 k 



m(T) = 



N 



Wo = -2(N 2 - l)V / ^ ln[l - exp(-£/?)] 
, dk k 2 ln[l - exp(-E/3)] 



{N 2 - 1)V 



2 3 / 2 3 7T 



g 2 (T) T T -> oo. (94) 



Thus, in the high-temperature limit, m(T) approaches a 
standard magnetic mass, to ~ g 2 (T)T. 



{N 2 - l)V/3 f dk (fc 4 - to 4 ) 



3tt 2 J E [exp(E(]) - 1] ' 
o 



(101) 
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where E = E(k) = (k 2 + ^-) 1//2 . We add the term W-2 where v = (x 2 + -^z) 1 / 2 ■ To leading order in r\ this gives 
and obtain to order g° the free energy per unit volume, 
w = W/V, 



2 /3 m 4 1 f dk (k 4 — m 4 ) 



2Ng 2 3tt 2 J E [exp(S/3) - 1 

(102) 

Here m = m(T) is the solution (|94| of the gap equation. 



X. EQUATION OF STATE AT HIGH 
TEMPERATURE 



We now evaluate w in the high-temperature limit, 
where g(T) is small, and our expansion should be reli- 
able. From (l94l we obtain 



w - 
where 



(» J -D(^0 + ^„))r3. ,103) 

oo 

dy (y 4 - v) 



K(n) = 



u (exp u — 1) ' 



and 



V 



u ^(y 2 + ^y/ 2 , 



m\ 4 (N g 2 {T)\± 



(104) 
(105) 

(106) 



. T / V 2 3 / 2 3 tt. 

is a small parameter. 

The lowest-order expression for K(rf) at small ry is ob- 
tained by setting r\ = 0, 



y 3 



exp y — 1 



7T 

15 : 



(107) 



which gives the Stefan-Boltzmann free energy. To evalu- 
ate the leading correction to K(0), we use 



dK 
dr] 



oo 

f d y (-r J --r 

J Vu (exp u — 1 



4 

y -v 



— 1) ' 2 u 3 y 2 (expit — 1) 
(y 4 - 77) exp u 



2 it 2 y 2 (exp u — l) 2 
We pose y = r] x l 4 x and obtain 

dK 
dr] 



(108) 



dx 



( L- 

\v fexpfji 1 / 4 ?;) 



[exp^ 1 / 4 ^) _ 1] 

x 4 - 1 
2 v 3 x 2 [exp^ 1 / 4 !,) - l] 

77 1 / 4 (x 4 — 1) exp(?7 :L / 4 w) 
~2v 2 x 2 [exp^ 1 / 4 ^) - ip 



(109) 



ftg -1 > / 1 x 4 — \ 

dr] rf-l 4 J 


-3 7T 



/ 1 3? 1 \ 

V t> 2 w 4 a; 2 / 



4-2V2.nl/4- 



(110) 



This identity, K = K(0) + Jq %, yields 



IT 4 7T T, 3 ' 4 



15 2V2 
We obtain for the free energy 



(111) 



(N 2 - 1) ( 



TV 3 CJ 6 7T 2 TV 3 £ 6 



2 7 3 3 7T 4 45 2 5 3 4 7T 4 
2 ^3 

n 6 m) t 3 

10,368 7T 4 



rpo 



g (T) J T , (112) 



which gives the leading correction to the Stefan- 
Boltzmann limit from the cut off at the Gribov horizon. 

The equation of state of the gluon plasma follows from 
the thermodynamic formulas for the energy per unit vol- 
ume and pressure, 



dw 
"dp' 



w 



(113) 



and entropy per unit volume, s = ^jr-. To calculate 



the energy density, we use — /?§§ 



T 
rpdg_ 
± dT 



/5-function 

= 0(g 3 ) 1 which is of higher order. We thus obtain for 
the energy density and pressure at high temperature, 



3p = (N 2 - 1) 



iV 3 



15 3,456 7T 4 



g\T)\ T , (114) 



with s = 

This provides the leading correction to the equation of 
state of the gluon plasma at high T that comes from the 
cut-off at the Gribov horizon. Numerically it is small, 
whereas the correction of order g 6 is divergent when 
calculated perturbatively [ll|. To this result must be 
added the perturbative contributions, including resuma- 
tions [43^ . Since standard, resummed perturbation the- 
ory diverges at order g 6 , which is precisely the order of 
the correction we have found, the result obtained here is 
consistent with standard perturbative calculations. 



XI. CONCLUSION 

We have developed a semi-perturbative method of cal- 
culation. It uses the systematic, perturbative expansion 
of local, renormalizable quantum field theory to calcu- 
late all quantities, including the gap equation. The gap 
equation is then solved non-perturbatively for the Gribov 
mass m(T). 
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We have have applied this method to calculate the 
leading correction to the Stefan-Boltzmann equation of 
state of the gluon plasma at high T that comes from the 
cut-off at the Gribov horizon. It is of order g 6 . How- 
ever we have not attempted to calculate the additional 
perturbative corrections to the Stefan-Boltzmann law to 
this order [43| . Significantly, g 6 (T) is precisely the order 
at which resummation of ordinary perturbation theory 
fails due to infrared divergences, as was first shown by 
Linde [ll|. Thus our calculation does not contradict the 
finite results from ordinary perturbation theory and its 
resummation, but instead gives a finite result precisely 
where ordinary perturbation theory breaks down. 
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APPENDIX A: ACTION IN TERMS OF 
HORIZON FUNCTION 

In this Appendix we integrate out auxiliary fields to ex- 
press the action in terms of the non-local "horizon func- 
tion" that effects the cut-off at the Gribov horizon. It is 
non-local in space, but local in time. 

We start with the action (Ti5|) . It contains cross terms 
that involve u> but none with to. These cross terms are 
cancelled by an appropriate shift of to. We then inte- 
grate out the quartet Wq 6 , Wg b , y)g b , (p^ b that appears in 
the action in the expression 

I = J d D x (^ c M Q Vo C " u>g c M ab u, b c ), (Al) 

where M ab = —diDf b is the Faddeev-Popov operator. 
This gives 

J dujoduiodipodfio exp(—Io) = 1, (A2) 

because from the fermi ghost pairs Wq 6 ,Wq 6 , we get 
(detAf)^ _1 and from the bose ghost pairs (/?g 6 ,<^g h we 
get l/(dctM) N _1 , which cancel. (This shows that the 
auxiliary quartet ipo , u)q , uiq , ipo plays no dynamical role 
but it was included because it may be useful to control 
Lorentz invariance.) 

We now integrate over the remaining auxiliary fermi 
ghost pairs and obtain 

J dufduf exp(wf , M ab wf) = (detM)^ 2 - 1 ^- 1 ). 

(A3) 



The remaining auxiliary bose ghost pairs appear only in 
the action 

= J d D x [<f¥M ah <p* + 7 1/2 9/*4(pf - ??*)], 

where we have integrated by parts. The integral over 
these ghosts is done by completing the square with the 
result 

/ dpfdyf expi-S^) (A5) 
= (detM)-( N2 -^ D -V exp(- 7 5 w ), 

where 

Shi = J d D x gf abc A b {M- l ) ad gf dec Al (A6) 

The powers of dct M from the fermi- and bose-ghost in- 
tegrations cancel, so the net result from integrating over 
all the auxiliary ghost pairs is simply exp(— 7 Shi)- 

The only remaining dependence on 7 occurs in the con- 
stant term in C m , CORRECT: eq. (gSJ). We combine S hl 
with the constant term and obtain the non-local horizon 
function 

S h =fd D x[ gf abc A b (M~ 1 ) ad gf dec A e t 

— (N 2 — 1)(D — 1) ]. (A7) 

The integal over the Faddeev-Popov ghosts, c and c, gives 
det M as usual, and the integral over the Nakanishi- 
Lautrup field gives Coulomb gauge condition, 8{diAi). 
The partition function is now coincides with (|27[) . More- 
over, the condition (|51[) that m, or 7, be a stationary 
point of the free energy W = In Z now reads 

dW 

= -(S h ) = 0, (A8) 

which coincides with ((28)) and (|30| . This establishes 
that the theory described by local, renormalizable La- 
grangian (|48j) supplemented by the gap equation (|5"T|) is 
equivalent to the standard Coulomb gauge with a cut-off 
at the Gribov horizon, as asserted in sect. II. 



APPENDIX B: CANONICAL COULOMB-GAUGE 
HAMILTONIAN 

In this Appendix we express the action in terms of 
the 3-dimensionally transverse field Af and its canonical 
momentum irf . 

Although Sh is non-local in space, it is local in time, 
because the 3-dimensional Faddeev-Popov operator M, 
and thus also its inverse M _1 , is local in time. Moreover 
Sh depends only on A4 but not on its time derivative Ai . 
The latter appears only in Sym = J d D x \F 2 V where 
F% = 2F& + F 2 and F Ql = Ai - D % A . This allows 
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us to introduce canonical momenta, which are the color- 
electric fields 7rf, by the Gaussian identity 

exp(-Sy M ) = / dm exp ( - J d D x [ i?if (Bl) 

- DiA ) + -Ffj } ), 

Integration over Aq imposes Gauss's law in the form 
d(DiTTi), which gives 

Z = J dAid-Ki det M 6(DiWi) SidiAi) (B2) 

x exp ( - J d D x [ ^7T? + m l A l + ^F? ] - 7 S h ) , 

We now separate transverse and longitudinal parts, A, = 
Af — did and m = nf — dill, where diAj = diirf = 0, 
and U is the color-Coulomb field, and obtain 

J dA.dm detM 8(D t m) <5(<Mi) 
= J dA T dn T dU detM 6 (Dim), (B3) 

where we have used f g?<t 8(dfa) = const. To solve 
Gauss's law we write 

Dim = D t (A T )(nf - diU) 

= gAf x n? + M(A T )U, (B4) 

where M(A T ) = —Di(A T )di is the 3-dimensional 
Faddeev-Popov operator. The Faddeev-Popov determi- 
nant gets absorbed by the identity 

det M J dU S(gAj x nj + MU) = 1. (B5) 

The last (^-function expresses Gauss's law, and fixes the 
color-Coulomb potential to its physical value, 

U = U phys = M- l p, (B6) 



where p = —gAf x irf is the color-charge density of the 
dynamical degrees of freedom. This gives the partition 
function in canonical form 

Z = J dA T diT T exp ( - J d D x [ ^tt? + mJAJ 

+ ±F*(A T )}- 1 S h (A T )), (B7) 

where m = ^7 ~ diU p h ys - This is a canonical sys- 
tem, and thus unitary, for all real 7, and stable for 
7 > 0. The canonical Coulomb-gauge hamiltonian of 
Christ and Lee [44] is modified here by the addition of 
the term jS h {A T ). 

The functional measure J dA T exp(-jSh), supple- 
mented by the horizon condition (|A8[) that fixes 7 = 
Iphys ■, were derived in [3l| as a representation of the func- 
tional integral restricted to the Gribov region, 

j dA T exp(-j phys S h ) = J dA T . (B8) 

u 

This yields the partition function in terms of the standard 
Coulomb-gauge first-order action 

Z = J dA T dTT T exp ( - J d D x [ + iwfAj 
n 

+ ±F*(A T )}), (B9) 

where m ; = n f ~ diUphysi and the integral over A T is 
restricted to the Gribov region 57, defined in sect. II. 

Whereas unitarity in Coulomb gauge holds in each or- 
der in the perturbative expansion, Lorentz invariance is 
violated at tree level. According to the equivalence with 
the operator formalism, Lorentz invariance is restored 
only when the non-perturbative gap equation is satis- 
fied, but it is not attempted to verify this by explicit 
calculation in the present article. 
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